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Abstract 

We study the peculiarities of stretching of globular polymer macromolecules 
in a disordered (crowded) environment, using the model of self-attracting self- 
avoiding walks on site-diluted percolative lattices in space dimensions d = 3. 
Applying the pruned-enriched Rosenbluth chain-growth method (PERM), 
we construct the phase diagram of collapsed-extended state coexistence when 
varying temperature and stretching force. The change in shape characteris- 
tics of globular polymers under stretching is analyzed as well. 
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1. Introduction 

Long flexible polymer macromolecules in a good solvent form crumpled 
coil configurations which are perfectly captured by the model of self-avoiding 
random walks (SAW) on a regular lattice [lj]. This regime holds at tempera- 
tures T well above the so-called G-point. When lowering the temperature, the 
effect of monomer-monomer attraction grows and the polymer radius shrinks. 
At T = Tq, a crossover occurs from high-temperature SAW behaviour to the 
G-statistics. At this particular temperature polymers in d = 3 dimensions 
behave effectively as simple random walks (RW). Below the Q-temperature, 
the entropic effects, which make the polymer chain swell, are overcome by 
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interaction energy and a collapse to the globule regime occurs. The coil- 
globule transition is considered to be of second order in the sense that 
the density of an infinite globule is zero at T = Tq and increases continuously 
when further lowering the temperature. 

The coil-globule transition is of interest in various respects, being deeply 
connected with problems like protein folding and DNA condensation. The 
properties of polymers in the vicinity of the 0-point can be successfully 
studied on the basis of self-attracting self-avoiding walks (SASAW), where a 
nearest-neighbor interaction is included: an attractive energy e between two 
neighbor sites is introduced. The coil-globule transition of flexible polymers 
has been so far the subject of numerous studies 0, H, 0, 0, @, 0]. Recent 
numerical estimates give T e (d = 3) = 3.717(3) (if. 

In studying the folding dynamics and transport properties of proteins, 
an important role is played by global shape properties of a typical polymer 
configuration. The asymmetry of polymer shape can be characterized, e.g., 
by the so-called averaged asphericity (Ad) jil, [1], which takes on a maximum 
value of one for a completely stretched configuration, and equals zero for a 



spherical form. It was realized experimentally [10J, that the majority of 
globular proteins are characterized by an asphericity value (Ad) ~ 0.1, thus 
being almost spherical. 

In polymer physics, of great importance is the understanding of the be- 
haviour of macromolecules in the presence of structural disorder. In partic- 
ular, related problems have been raised recently in studies of protein folding 



in the natural cellular environment [12|, [13[ . Real biological cells can be de 



scribed as a very crowded environment built of the biochemical species, which 
occupy a large fraction of the total volume. In the language of lattice models, 
the crowded environment with structural obstacles can be considered as a dis- 
ordered lattice, where some amount of randomly chosen sites contains defects. 
Of particular interest is the case, when the concentration p of lattice sites 
allowed for the SAWs equals the critical concentration p c (d=3) = 0.31160 



14j and the lattice becomes percolative. It is established that the value of 



the 6-temperature is lowered due to the presence of disorder [15|, [16|, llTJ, [18 
numerical estimates give Tq(cL = 3) = 0.71(2) 

The recent progress in experimental techniques makes it possible to moni- 
tor the behaviour of various polymers under tension and stress. In particular, 
appl ying a force on an isolated protein, the unfolding of the giant titine pro- 



tein [19] and stretching of collapsed DNA molecules 20] have been studied. 



Of special interest in biophysics is the stretching of globular polymers below 
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the 0-point. Force not only influences the structural properties of polymers, 
but also may introduce a new completely stretched state which is otherwise 
not accessible. The properties of force-induced transitions in polymers have 
been studied intensively 2l|, |22|, |23_|, 125] ■ The response of a polymer in 



crowded media to the stretching force modelled by the SASAW model on a 



percolative lattice has been analyzed recently in [26| and [18fl . The interesting 



question about how the shape properties of almost spherical polymer globules 
are influenced by stretching remains, however, completely unresolved. 

The aim of the present study is to apply numerical simulations to analyze 
the properties of SASAWs on site-diluted lattices at the percolation thresh- 
old under applied external stretching force in space dimensions d = 3. We 
analyze the effect of applied force on the phase transitions between collapsed, 
extended and stretched phases and estimate the influence of stretching on 
the shape parameters of globular proteins in crowded environments. 



2. The method 

We consider site percolation on regular lattices of edge lengths up to 
L max =200 in d — 3. Each site of the lattice was assigned to be occupied with 
probability p c and empty otherwise. To obtain the backbone of a percolation 
cluster on a given disordered lattice, we apply an algorithm explained in 
detail in our previous papers [27 . 

To study SASAWs on the backbone of percolation clusters, we apply 
the pruned-enriched Rosenbluth method (PERM) taking into account 
that the SASAW can have its steps only on the sites belonging to the back- 
bone of the percolation cluster. PERM is based on the original Rosenbluth- 
Rosenbluth (RR) method 28( and enrichment strategies [29]. The polymer 
grows step by step, i.e., the nth monomer is placed at a randomly chosen 
empty neighbor site of the last placed (n — l)th monomer (n < N, where N 
is the total length of the chain). The growth is stopped, if the total length 
of the chain is reached. A weight 



W n ][ 



mie 



(1) 



1=2 



is given to each sample configuration at the nth step, where mi is the number 
of free lattice sites to place the Ith monomer, Ei denotes the energy of the 
/-step chain (E\ = zi-e with e being an attractive energy between two nearest 
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neighbors and z\ the number of nearest neighbors contacts for a given chain) 
and ks is the Boltzmann constant. In what follows, we will assume units in 
which e = —1, ks = 1. 

The configurational averaging for any quantity of interest then has the 
form: 

(•••> = ^E^ onf ---> ^v = E^ onf - ( 2 ) 

Zjjy — — 
conf conf 

This method is particularly useful for studying G-polymers, since the Rosen- 
bluth weights of the statistically relevant chains approximately cancel against 
their Boltzmann probability. 

Population control in PERM suggests pruning configurations with too 
small weights, and enriching the sample with copies of high-weight configu- 
rations [5] . These copies are made while the chain is growing, and continue to 
grow independently of each other. Pruning and enrichment are performed by 
choosing thresholds W< and W> depending on the estimate of the partition 
sum of the n-monomer chain. If the current weight W n of an n-monomer 
chain is less than W^, a random number r=0, 1 is chosen; if r=0, the chain 
is discarded, otherwise it is kept and its weight is doubled. Thus, low-weight 
chains are pruned with probability 1/2. If W n exceeds W>, the configuration 
is doubled and the weight of each copy is taken as half the original weight. 

In stretching studies, one end of the chain is subjected to an external 
force F acting in a chosen direction, say x, while the other end (the starting 
point) is kept fixed. The stretching energy E s arising due to the applied force 
for an n-step trajectory is given by: E s = —F ■ x, where x = x n — x denotes 
the x-component of the distance from the starting point. The Rosenbluth 
weight factor W n is here taken to be: 

1=2 

For estimations of quantities of interest we have to perform two types 
of averaging: the first over all polymer configurations on a single percola- 
tion cluster according to (j2J); the second average is carried out over different 
realizations of disorder, i.e., over all percolation clusters constructed: 

i M 

£<•••>, (4) 

1=1 
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Figure 1: Specific heat per monomer of a SASAW with N = 90 steps in d = 3 under 
stretching force F as a function of temperature. Left: pure lattice, right: backbone of 
percolation cluster. Squares: F = 0.2, diamonds: F — 0.4, triangles: F = 0.6. 



where M is the number of different clusters and the subscript i means that 
a given quantity is calculated on the cluster i. 

3. Results 

The properties of systems in the vicinity of a second-order phase transi- 
tion can be studied by analyzing the peak structure of the specific heat Cy as 
a function of temperature indicating crossovers between physically different 
states. In the case of a polymer system, this corresponds to the transition be- 
tween globule and coil regimes. Cy can be expressed via energy fluctuations 
as follows: 

Cv{T) = 1 — 2 {W)-W)- (5) 

To study the 0-transition of SASAWs, when the external stretching force 
is acting in the environment, we are working in the "constant-force" ensem- 
ble. Fixing the value of F, we study the specific-heat behaviour (Fig. [I]). 
Analyzing the peak structure of the specific heat, we immediately conclude, 
that increasing the value of F leads to decreasing the transition tempera- 
ture. With increasing F, the averaged energy of the chain decreases, because 
applied force stretches the polymer globule. The value of the transition tem- 
perature is thus shifted by the presence of force. 

For finite chain length N, the temperature defined by the position of 
the specific- heat maximum T^ X (N) is well below the collapse transition 0- 
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Figure 2: Phase diagrams of a SASAW under applied force F in c£ = 3. Left: pure lattice 
right: backbone of percolation cluster. The filled triangle presents the result of Ref. 23] 
F = 1.5, Te ~ 2.46. 



temperature. This finite-size deviation of T™ X (N) from T® obeys scaling 
behaviour with N: 

T™ x (A0-T e ~a-A^ + A (6) 

where a, b are constants and ^e(ci = 3) = 1/2 is the size exponent of a SAW 
at the 0-point. Our estimates for Tg c in the presence of a stretching force are 
obtained by least-square fitting of (JB} in Ref. [l8[. Results are presented in 
Fig. |2] in the form of a phase diagram of transitions from globular to extended 
states. 

The measure of the shape properties of a specified spatial conformation 
of a polymer chain can be characterized [9| in terms of the gyration tensor 
Q with components: 



Qi 



l N 



x 



CM 



)(4 



r 



cm '■ 



), i,j = l,...,d, 



n=l 



where x 3 n is the jth coordinate of the position vector of the nth monomer of 
a polymer chain (n = 1, . . . , N), and x J CM = J2n=i x rJN is the coordinate 
of the center-of-mass position vector. The extent of asphericity of a polymer 
configuration can be characterized by the quantity defined as Jsf] : 



(A,, - A) s 



d(d 



1 X- i/H - 



d TrQ 2 
d - 1 (Tr Q) 



2 : 



(7) 
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Figure 3: Averaged asphericity of an N — 90-step SASAW as function of applied force F 
in d — 3. Left: pure lattice (T = 1.8), right: backbone of percolation cluster (T = 0.2). 
Lines are guides to the eyes. 



with \i being the eigenvalues of the gyration tensor, A = TrQ/cf, and Q = 
Q — A I (here I is the unity matrix). This universal quantity equals zero for 
a spherical configuration, where all the eigenvalues are equal, A, = A, and 
takes a maximum value of one in the case of a rod-like configuration, where 
all the eigenvalues equal zero except of one. Thus, the inequality < A d < 1 
holds. 

In Fig. E]we present the averaged asphericity of SASAWs, giving informa- 
tion about the internal structure of the polymer configuration under applied 
force, at a temperature well below the 0-point for the cases of a pure lattice 
and the backbone of a percolation cluster. Note that in the absence of force, 
(A 3 ) for SAWs on a pure lattice is very close to zero, whereas in the disordered 
case, due to the complicated structure of the underlying percolative lattice, 
globular configurations are more elongated with larger (A^ c ) values. At small 
F, a polymer chain is still in the compact folded state and is just slightly 
oriented along the force direction. Under inreasing forces, the polymer chain 
takes on a conformation similar to the extended (swollen) structure. Note, 
that completely stretched states, corresponding to (Ad) — 1 can be obtained 
only in the pure case and are not accessible on the percolative lattices due 
to the complicated fractal structure on the underlying percolation cluster. 
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4. Conclusions 

We studied self-attracting self-avoiding walks on disordered lattices in 
space dimensions d — 3, modelling flexible polymer macromolecules in porous 
environment. We considered the special case, when the concentration of dis- 
order is exactly at the percolation threshold, so that an incipient percolation 
cluster of sites, allowed for SAWs, emerges on the lattice. Keeping one end 
of a SASAW trajectory on the backbone of a percolation cluster fixed, we ap- 
plied a stretching force F, acting in some chosen direction (say, x). Based on 
our numerical simulation data, we constructed phase diagrams of collapsed 
and extended states coexistence. The behaviour of averaged asphericity of 
globular SASAWs on a percolation cluster under a stretching force was ana- 
lyzed as well. As expected, in the case of a pure lattice, stretching leads to 
increasing the asymmetry indicating at some critical value of F a transition 
to a completely stretched state with (Ad) ~ 1. For SASAWs on percolative 
lattices, the completely stretched states are not accessible due to the compli- 
cated fractal structure of the underlying percolation cluster even under very 
strong stretching. 
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